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ABSTRACT 

Context. Magneto-rotational instability (MRI) has been suggested to lead to a rapid growth of the magnetic field in core collapse 
supernovae and produce departures from spherical symmetry that are important in determining the explosion mechanism. 
Aims. We address the problem of stability in differentially rotating magnetized proto-neutron stars at the beginning of their evolution. 
Methods. To do this we consider a linear stability taking into account non-linear effects of the magnetic field and strong gravity. 
Results. Criteria for MRI are derived without simplifying assumptions about a weak magnetic field. In proto-neutron stars, these 
criteria differ qualitatively from the standard condition dfl/ds < where £1 is the angular velocity and s the cylindrical radius. If 
the magnetic field is strong, the MRI can occur only in the neighbourhood of the regions where the spherical radial component of the 
magnetic field vanishes. The growth rate of the MRI is relatively low except for perturbations with very small scales which usually are 
not detected in numerical simulations. We find that MRI in proto-neutron stars grows more slowly than the double diffusive instability 
analogous the Goldreich-Schubert-Fricke instability in ordinary stars. 
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1. Introduction 

There is a growing amount of evidence that core-collapse su- 
pernovae are asymmetric and that the core-collapse mechanism 
itself is responsible for the asymmetry (see Buras et al. 2003, 
Akiyama et al. 2003 for more details). Several possibilities are 
explored to account for this observed asymmetry. One is associ- 
ated with the influence of rotation on convection, which seems 
to be inevitable during the early evolution of proto-neutron stars 
(PNS). Convective motions in PNSs are very fast (~ 10 8 — 10 9 
cm/sec) and, therefore, the convective turnover time is short, 
~ 1 — 10 ms (see, e.g., Burrows & Lattimer 1986). Nevertheless, 
if angular momentum is conserved, the collapsing core can spin 
up to very short periods ~ 5 — 10 ms and generate strong differ- 
ential rotation. Such fast rotation modifies convection and makes 
convective motions anisotropic and constrained to the polar re- 
gions (Fryer & Heger 2000, Miralles et al. 2004). This mecha- 
nism is a natural way to create anisotropic energy and momen- 
tum transport by convective motions, that only requires that the 
angular velocity be of the order of the Brunt-Vaisala frequency. 

The other possibility to create asymmetry is the effect of 
jets (see, e.g., Khokhlov et al. 1999, Wheeler et al. 2002). Even 
though the mechanism of jet formation is still unclear, it seems 
that MHD jets are common in systems where a central body ac- 
cretes matter with angular momentum and magnetic field (see, 
e.g., Meier at al. 2001), and a core-collapse supernova is such 
a system. Calculations have established that nonrelativistic axial 
jets originating within the collapsed core can initiate a bipolar 
asymmetric supernova explosion that is consistent with observa- 
tions (Hwang et al. 2000). 

Another way to generate asymmetry is associated with the 
magnetic field that can be an important ingredient of the ex- 
plosion mechanism (Bisnovatyi-Kogan 1971, Kundt 1976). The 
toroidal magnetic field can be amplified by differential rota- 



tion to such high values that it becomes dynamically important 
(Ardelyan et al. 2005). The effect of the magnetic field on asym- 
metry of supernovae was considered by Wheeler et al.(2000, 
2002) who found that it is possible to produce both a strong 
toroidal field and an axial jet. Two-dimensional MHD simula- 
tions of core collapse indicate that the shape of shock waves and 
the neutrinosphere can be modified by the effect of the magnetic 
field (Kotake et al. 2004, Takiwaki et al. 2004). 

The possible presence of a magnetic field and differential 
rotation in a core-collapse supernova favours magnetorotational 
instability (MRI), which can enhance turbulent transport and 
amplify the magnetic field. This instability was considered by 
Akiyama et al.(2003) in the context of core collapse. The au- 
thors argued that instability must occur in core collapse and 
that it has the capacity to produce fields that are sufficiently 
strong to affect, if not cause, the explosion. Thompson et al. 
(2005) constructed one-dimensional models, including rotation 
and magnetic fields, to study the mechanism of energy deposi- 
tion. They explored several mechanism for viscosity and argue 
that turbulent viscosity caused by the MRI can be most effec- 
tive. Numerical simulations provide contradictory conclusions 
regarding the importance of MRI in core collapse. Moiseenko 
et al. (2006) claim that MRI has been found in their 2D simu- 
lations, and that it is responsible for a strong amplification of 
the poloidal magnetic flux. However, what these authors call 
MRI is different to standard MRI considered by Velikhov (1959) 
(see also Balbus & Hawley (1991)). For instance, the instability 
found by Moiseenko et al (2006) starts to develop only when the 
ratio of the toroidal and poloidal fields reaches a value of ~ a few 
tens. On the other hand, the onset of standard MRI in 2D does 
not depend on the toroidal field at all. The dependence on the 
ratio of the toroidal and poloidal field is more typical for Tayler 
instability (Tayler 1973), which is more relevant to the topol- 
ogy of the magnetic field than differential rotation. Therefore, 



2 



V. Urpin: Magnetorotational instability in proto-neutron stars 



it is possible that Moiseenko et al. (2006) incorrectly identify 
instability, and MRI does not occur in their simulations. Apart 
from that, Moiseenko et al. (2006) attributed a rapid growth of 
the toroidal and poloidal fields to a dynamo driven by the mag- 
netorotational instability. This also rise some doubts because of 
Cowlings's anti-dynamo theorem which states that an axisym- 
metric dynamo cannot exist (see, e.g., Shercliff 1965). Two- 
dimensional simulations of core collapse with a strong magnetic 
field have been performed by Sawai et al. (2005, 2008). They 
found that the magnetic field can play an important role in the 
dynamics of the core only if the poloidal field of the progeni- 
tor is strong enough (<~ 10 12 — 10 13 G), but MRI was not seen 
in the considered models. The magnetic field is amplified mainly 
by field compression and field wrapping in these simulations. On 
the contrary, Shibata et al. (2006) claim that they found MRI in 
their simulations of magnetorotational core collapse in general 
relativity. These authors paid attention to resolution in order to 
resolve unstable MRI modes and they claim that amplification of 
the magnetic field in the considered models is caused by MRI. 
Note, however, that the poloidal field obtained in their models is 
of the order of that estimated from conservation of magnetic flux 
in core collapse, and a more refined analysis is required to de- 
termine the mechanism of amplification. Fryer & Warren (2004) 
argued that it is difficult to produce magnetic fields in excess 
of 10 14 G even if MRI occurs in core collapse because rota- 
tion is not sufficiently fast. A detailed study of the magnetoro- 
tational core collapse has been performed by Obergaulinger et 
al. (2006a, 2006b, 2009). The initial magnetic field was purely 
poloidal in their models with a strength ranging from 10 10 to 
10 13 G. Such fields are much higher than those estimated to ex- 
ist in realistic stellar cores, but the authors wanted to investigate 
the principal effects of a magnetic field. The initial magnetic 
field is amplified by differential rotation in these simulations, 
giving rise to a strong toroidal component. The poloidal compo- 
nent grows mainly by compression during collapse and does not 
change significantly after core bounce. The authors also found 
that extended regions exist where the criterion of MRI is satis- 
fied at various epochs. However, the growth rate of this instabil- 
ity is typically too small, except for a few models with a strong 
initial magnetic field B = 10 12 G. 

In this paper, we study the effect of MRI on core-collapse 
supernova. Since the magnetic field can be sufficiently strong, 
the criterion of instability is derived taking into account terms 
depending on the Alfven frequency. We derive the criterion that 
applies to any rotation profile but a special consideration is made 
for the case of shellular rotation that often is used to mimic rota- 
tion of proto-neutron stars. We address only the axisymmetric 
instability because numerical simulations of a magnetic core- 
collapse are usually done in 2D. The main goal of this study 
is to show that the effect of MRI on proto-neutron stars often is 
overestimated. 



2. The growth rate of convective and 
magnetorotational instabilities 

We assume that the initial PNS is restricted by the radius of neu- 
trino sphere. The PNS has a high-entropy mantle, so that the 
outer part of the star is initially at a relatively large radius. It 
takes a few tenths of a second for the neutrinos in the high- 
entropy mantle to leak from the star and for the mantle to col- 
lapse to the canonical radius. In this paper, we study MRI in rel- 
atively deep layers of the PNS where the density is comparable 
to (or higher than) the nuclear density. 



Consider a PNS rotating with angular velocity SI = Sl(s, z); 
(s, (p, z) are cylindrical coordinates. We explore the Boussinesq 
approximation and assume that the magnetic energy is small 
compared to the thermal one. In the unperturbed state, the star is 
in hydrostatic equilibrium, 

— = G+—(VxB)xB, G = g + n 2 s (1) 
p 4np 

where g is the gravity. Generally, the magnetic field B has both 
toroidal B v and poloidal B p components. If the field has a com- 
ponent parallel to Vfi, the azimuthal field in the unperturbed 
state increases with time by winding up the polodal field lines. 
If the magnetic Reynolds number is large, then 

B v {t) = B v (0) + s(B p -Vn)t. 

We consider MRI assuming that the basic state is quasi- 
stationary. This is justified if the growth rate of the MRI (~ J!) 
is greater than the inverse time-scale on which the basic state 
evolves. This inverse time-scale can be estimated as \p/p\ ~ 
s\\7VL\(pm /p) ~ (pm/p)Q where pu is the magnetic pres- 
sure. Therefore, the condition of quasi-stationarity is satisfied 
if Pm Ip < 1. that is fulfilled in core collapse supernovae. 

We consider axisymmetric short-wavelength perturbations 
with spatial and temporal dependence exp(7i — ik ■ r) where 
k = (k s ,0, k z ) is the wave-vector. Small perturbations will be 
indicated by a subscript 1 . Then, the linearized MHD equations 



read in a short-wavelength approximation 

7«i + 2fl x V! + e^sfuiVO) = — ^ + G — + 

P P 

-^—[k(B ■ £?i) - Bi{k ■ £?)] , (2) 

fe • vi = , (3) 
7Bi = -i«i(fc-B) + ae v ,(Bi-Vn), (4) 
k ■ B x = , (5) 



where e v is the unit vector in the ^-direction. We neglect kine- 
matic and magnetic viscosities since they are small in PNSs. 
The term oc V x B in Eq. (2)) is also neglected. In the short- 
wavelength approximation, the latter is justified if B v satisfies 
the inequality B v < (L/X)B p , where L is the length scale of 
unperturbed quantities and A = 2n/k. 

The matter is assumed to be in chemical equilibrium, thus 
the density is a function of the pressure p, temperature T and lep- 
ton fraction Y. In the Boussinesq approximation, perturbations 
of the pressure are small and, therefore, p\ can be expressed in 
terms of T\ and Y\, 

pi « -p/?(Ti/T) - pSY u (6) 

where j3 and S are the coefficients of thermal and chemical 
expansion; (3 = — (<91np/<91nT) p y, S = — (din p/dY) p x. 
Since diffusive timescales are long compared to dynamical ones 
(Miralles et al. 2000), the linearized transport equations read 

Ti - vi ■ AVT = 0, (7) 
Y + Vl ■ VY = 0, (8) 

where AVT = [dT/dp) s ^P - VT is the super-adiabatic 
temperature gradient. 
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The dispersion equation corresponding to Eqs. (2)-(8) is 
7 4 + 6 27 2 + 6 = 0, (9) 
where 

b 2 = 2u? A + oj 2 + q 2 , b = uj\{uj 2 a + u 2 c + q 2 - 40 2 fc 2 /fc 2 ). 
Here ui A = (k ■ B)/ ^Airp is the Alfven frequency and 

q 2 = (k 2 z fl 2 e - sk s k z fl 2 z )/k 2 , u 2 c = -C-[G- k(k ■ G)/k 2 }, 
where C = -((3/T)AVT+8VY, Q e is the epicyclic frequency, 

n 2 e = d(s 4 n 2 )/s 3 ds, and n 2 z = dn 2 /dz. 

Eq. (9) has an unstable solution if either b 2 or b is negative, 

or 

2u 2 A + lo 2 + q 2 < 0, uj\ + uj 2 + q 2 - m 2 {k 2 z /k 2 ) < 0. (10) 
In the limit of a weak field, these conditions yield 

LO 2 + q 2 <0, LO 2 + q 2 -4n 2 (k 2 /k 2 )<0. (11) 

The first inequality is the criterion of convection modified by ro- 
tation, and the second condition represents the criterion of MRI. 
The solution of Eq. (9) is 



2 2 1 / 2 , 2\ 1 

7 =-W4--K+<? )± 



iK 2 + , 2 ) 2 +4^ 2 | 



1/2 



.(12) 



If uja is small compared to the other frequencies, we have 



-U + q% 



73 2 4 



2 

-UJ A 



40 



21,2 



k 2 (co 2 +q 2 ) 



■ (13) 



The solutions 71,2 correspond to buoyant modes that cause con- 
vection and are unstable if ui 2 + q 2 < 0. The solutions 73,4 
describe the magnetorotational modes which can be unstable if 
the second condition (12) is satisfied. 

For the purpose of illustration, we plot in Fig. 1 the depen- 
dence of 7 2 on oja/Q for the unstable magnetorotational mode. 
Its growth rate is given by Eq.(12) with the upper sign. Even if 
stratification is negligible, the growth rate of this mode is typ- 
ically low and 7 ~ u>a in a weak magnetic field (u)a < O). 
Stratification can substantially decrease the growth rate and this 
is seen very well from the figure. For example, the growth 
rate decreases approximately by a factor of two if uj 2 /Vl 2 . At 
lo 2 /VL 2 w 1, stable stratification completely suppresses the mag- 
netorotational instability of the considered perturbations. 

In the convective zone where uj 2 + q 2 < 0, the quantity 
lo 2 + q 2 — 4D?(k 2 /k 2 ) is also negative and, hence, 734 < 
0. Therefore, the magnetorotational instability does not occur 
in convectively unstable regions, and these two instabilities 
are spatially separated if the magnetic field is weak (see also 
Obergaullinger et al. 2009). The value of io c is of the order 
of 1-10 ms _1 in collapsing cores (see Thompson et al. 2005). 
Therefore, rotation has an important impact on convection only 
if ft is of the same order of magnitude, ft s=a 1000 rad/s (Miralles 
et al. 2004). This value can be reached in PNS if rotation of the 
progenitor was very fast (Villain et al. 2004). On the contrary, 
MRI can have an important influence even if rotation is slower. 



(M 

a 




0.4 0.6 0.8 

a) A / Q. 

Fig. 1. The dependence of 7 2 onwA/OforO oc l/sandfc z /fc = 
V2/2 and for three values of the ratio w 2 /0 2 = 0.0, 0.5, and 
0.8. 

3. The criteria of magnetorotational instability 

The condition of MRI (second inequality (10)) depends on the 
direction of k and can be written as follows 



bo 



A 



fc 2 



D 



k s k z 



k 2 



(14) 



where 

A = sVL 2 + cos 2 au\o + C Z G Z , E = sin 2 aLu\ + C S G Z , 
D = sn 2 z - sin 2a uj 2 A0 - C S G Z - C Z G S . 

In these expressions, we denote 

k 2 B 2 



ujo 



2 -C G 



2 



Airp 



dn 2 



(15) 



a is the angle between the magnetic field and the rotational axis, 
cosa = B z /^B 2 +Bl 

Since the dependence of F on the direction of k is simple, 
we can obtain that F reaches its minimum at 



k 2 



1± 



{A-Ef 



(A - Ef + D 2 



(16) 



The value of F corresponding to these k 2 / k 2 yields the follow- 
ing condition of instability 



sn 2 s +lo 2 + uj 2 aq ± y/D 2 + (A - E) 2 < 0. 



(17) 



By taking the curl of Eq. (1), it can be readily obtained that the 
condition of hydrostatic equilibrium leads to 



where 
G B 



[C x Gb — L] = C z Gbs — C s Gbz — L v 



(18) 



G 



Airp 



(V xB)xB , L = Vx(G B -G). 



We consider only the magnetic field satisfying the condition 
L <C sD, 2 that is approximately equivalent to the requirement 
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that the magnetic energy is small compared to the rotational en- 
ergy. Then, Eq. (17) can be further simplified to obtain 

sill + a, 2 +^ 2 ±{ ( s0 2 + col + u 2 AQ ) 2 + 
G z [s(C z n 2 s - C s n 2 z ) + u 2 A0 {C z cos 2a + C s sin 2a)] - 

sin 2 auj 2 A0 (Lu 2 + sQ 2 s )} < 0. (19) 

The two conditions of instability follow straightforwardly from 
the above expression: 

stf s + to 2 + uu 2 m < , (20) 
sG z (C z n 2 s -C s n 2 z )+ uj 2 aq [C z G z cos 2 a + C s G z sin 2a 

+ (C S G S -sn 2 )sm 2 a] >0. (21) 

Conditions (20) and (21) look like the Solberg-H0iland condi- 
tions (Tassoul 2000), but with additional terms due to the chemi- 
cal composition gradients and the magnetic field. If the magnetic 
field is weak and V7 = then Eqs. (20) and (21) yield 

s n 2 s + luI < o ; sG z (c z n 2 s - c s n 2 z ) > o . (22) 

The first criterion is similar to the Schwarzschild criterion for 
convection modified by rotation. However, the Schwarzschild 
criterion at B = involves the angular momentum gradient, 
whereas the criterion of magnetorotational instability depends 
on the angular velocity gradient, as was noted by Balbus (1995). 
If rotation is cylindrical with il 2 = 0, the second criterion 
(22) yields the standard condition of the MRI, il 2 < 0, since 
C Z G Z < in a convectively stable region. 

Taking into account that C z tt 2 - C s Vl 2 z = (C x Vtt 2 ) v , 
Eq.(21) can be transformed into 

sG z | C 1 1 Vfl 2 1 sin ip + lo ao [C z G z cos 2 a + C s G z sin 2a 

+ (C S G S - sn 2 s ) sin 2 a] > , (23) 

where ip is the angle between vectors C and Vf2 2 . If the Alfven 
frequency is small, then the criterion of MRI reads 

> V > -7T, (24) 

and is different from the usually used condition 51 2 < 0. 
Criterion (24) implies that the component of Vfi 2 perpendic- 
ular to C must have a negative projection on e s for instability 
whereas the component along C plays no role. In the general 
case when loao is comparable to (or greater than) other frequen- 
cies (generally, that is possible for very short wavelengths even 
if B is weak), criterion (23) can be complicated. For instance, if 
lj ao ^> O, then criterion (23) yields 

C Z G Z cos 2 a + C S G Z sin 2a + (C S G S - sft 2 ) sin 2 a > . (25) 

This condition depends on the direction of B and is different 
from the standard condition Vl s < 0. Note, for example, that 
this condition cannot be satisfied if B s = because in this case 
sin a = and C Z G Z < in a convectively stable region. 

4. Instability in core-collapse supernovae 

The occurence of MRI is sensitive to the rotation profile. It fol- 
lows from both theoretical modelling and analytic consideration 
that core collapse of a rotating progenitor results in differen- 
tial rotation of a protoneutron star (Zwerger & Miiller 1997, 
Dimmelmeier et al. 2002, Miiller et al. 2004, Obergaulinger et 
al. 2006a,b). Many studies of MRI model rotation of the col- 
lapsing core by a shellular profile with Q = O(r) where r is 



the spherical radius (see, e.g., Akiyama et al. 2003, Thompson 
et al.2005, Sawai et al. 2005). Such rotation can be justified 
if the progenitor rotates with an angular velocity that depends 
on r alone (Monchmeyer & Miiller 1989). Then, if the angu- 
lar momentum is conserved, it can be shown that rotation of 
the collapsing core is shellular at the beginning of evolution, at 
least (see, e.g., Akiyama et al. 2003). However, there are studies 
that assume a cylindrical rotation for the initial profile (see, e.g., 
Obergaulinger et al. 2006a,b). 

We assume that the angular velocity is low compared to the 
Keplerian one and little departure occurs from a spherical geom- 
etry. Then, C and G are approximately radial, and we have 

G s = -gsm6, G z = -gcos6, C S = C sin 6, C Z =C cos 8, (26) 

where 9 is the polar angle. The quantity C is negative in a con- 
vectively unstable region and positive in a stable region. 

Substituting expressions (26) into Eqs. (20)-(21), we obtain 

stt 2 + lj 2 + lj 2 A0 < ; (27) 
sloq cos 6*(cos #£! 2 — sin Oft 2 ,) + 
to A0 [lu 2 (cos 6 cos a + sin sin a) 2 + sfi 2 sin 2 a] < 0. (28) 

Condition (28) can be satisfied only if the angular velocity is 
higher than the Brunt-Vaisala frequency which is rather high 
in core-collapse supernovae. For example, according to calcu- 
lations by Thompson et al. (2005), the value of ujo decreases 
from <~ 10 3 to 10 2 s _1 as the radius increases from 50 to 200 
km in the gain region. The angular velocity can be comparable 
only if the progenitor rotated very rapidly with a period ~ 2 s. 
Such rotation can be achieved only in the most rapidly rotating 
stars, and our study does not address such stars. We concentrate 
on condition (28) that can be satisfied for much slower rotation. 
If rotation is approximately shellular as is often assumed and 
w Cl(r), then we have n 2 s — f2 2 sin# and Vt 2 = fl 2 cos9 
where fi 2 = dti 2 /dr. Then, Eq.(28) can be transformed into 

lu 2 cos 2 (0 -a)+r sin 2 9 sin 2 a ft 2 < 0. (29) 

It turns out that stratification strongly suppresses MRI for such 
a rotation profile because lu 2 > in convectively stable regions 
and, typically, uj 2 > Q 2 . 

In Fig. 2, we plot the critical value of rVL 2 that discriminates 
between stable and unstable regions as a function of 9 for several 
values of the angle a. The region above the corresponding curve 
is magnetorotationally unstable. For any given a and negative 
f2 2 , there exists a range of 9 where the instability can arise. If 
rf2 2 is greater than (or comparable to) to 2 then the instability 
occurs over a rather wide range of 9. If r|fi 2 | <C oj 2 (which is 
more typical for PNSs) then the instability arises only in a very 
narrow range of 9. This dependence can be easily understood 
from Eq. (29). Even if il 2 < 0, MRI occurs only in those regions 
where the first term on the l.h.s. of Eq. (29) is small. This occurs 
in the neighbourhood of the line 

9 = a± tt/2. (30) 

Condition (30) implies that the field line is perpendicular to the 
radius and, hence, at any magnetic topology, the instability oc- 
curs only near the region where B r = 0. How extended this 
region is depends on the relation between uj 2 and rVL 2 . 

There is no generally accepted point of view regarding topol- 
ogy of the magnetic field in core-collapse supernovae. For illus- 
tration, we consider the simplest configuration with 

B r = /(r)(l - 3 cos 2 (9), B e = F(r) sin0cos6>, (31) 
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Fig. 2. The dependence of the critical value of r|Sl 2 | that deter- 
mines the instability on the polar angle for a = — tt/3 (curve 
1). — tt/6 (2), 7r/6 (3), and tt/3 (4). The regions above the lines 
correspond to instability. 



where / and F are functions of the spherical radius that satisfy 
the divergence-free condition F = d(r 2 f)/rdr. This field can 
be matched to the dipole component outside the star (see, e.g., 
Urpin et al. 1994). Generally, the field can have a more complex 
geometry but a simple model (31) allows us to show qualitatively 
how the magnetic topology infuences the region of MRI. We 
have for the magnetic field (31) 



sm a 



cos a 



e s B 


sin 9 


Bp 


B 


e z B 


cos 9 


Bp 


B 



[J {I- 3 cos 2 9) + Fcos 2 9], 
[/(l - 3 cos 2 (9) -F sin 2 19]. (32) 



Substituting these expressions into Eq. (29), we obtain 

cj 2 (l-3cos 2 6») 2 +rf7 2 sin 4 6»[l-(l-x)cos 2 6»] 2 < 0, (33) 

where x = din //din r. Obviously, this inequality can be sat- 
isfied only if the angular velocity decreases with r. If SI 2 < 0, 
then criterion (33) is always fulfilled at cos 9 w l/y/3 (or at 
9 ~ 9q — 54°) except for the case x = —2. Hence, some region 
in the neighbourhood of this polar angle can also be unstable, 
and the opening angle of the unstable region depends on the ra- 
tio |rSl 2 /uJq I . For example, if this ratio is small then MRI occurs 
within the cone 9 + A9 > 9 > 9 - A9, where 



A9. 



3\/2 



LJ 



(34) 



5. Conclusion 



Interest in MRI in core collapse supernova is due to the fact that 
it can generate a strong magnetic field and produce significant 
departures from spherical symmetry. These departures are cru- 
cial for the explosion mechanism. In the present paper, we have 
derived criteria of MRI in proto-neutron stars taking into account 
the effect of a non-axial magnetic field. Criteria have been ob- 
tained in a form analogous to the Solberg-H0iland criteria but 
including terms containing the magnetic field. It turns out that 
the criterion of MRI in proto-neutron stars can differ from the 
standard condition dto/ds < even in a weak magnetic field 



due to strong gravity and the gradient of the lepton fraction. If the 
Alfven frequency ujao is small compared to the angular velocity 
then the instability occurs in the region where the component of 
Vfi 2 perpendicular to C = -(/J/T)AVT + SVY has a nega- 
tive projection on e s whereas the component VS1 2 parallel to C 
is unimportant for instability. In the case of slow rotation, when 
departures from sphericity are small (g 3> sSl 2 ), the criterion re- 
duces to a simple inequality dVl/d9 < (see, e.g., Urpin 1996). 
For instance, shellular rotation with SI = to(r) which is often 
used in modelling proto-neutron stars does not satisfy this con- 
dition. Therefore, MRI does not occur if the proto-neutron star 
rotates shellularly and the wavelength of perturbations is such 
that the Alfven frequency is smaller than SI. Only detailed calcu- 
lations of rotational core collapse can give the answer to whether 
the condition of MRI is fulfilled in proto-neutron stars and, gen- 
erally, this answer should depend on rotation of the progenitor. 
Note that the velocity of unstable perturbations is approximatelly 
perpendicular to the radius because gravity strongly suppresses 
motion in the radial direction and, likely, the radial turbulent 
transport should be suppressed when the instability saturates. 

However, MRI can arise in proto-neutron stars even if the 
neccesary condition that determines the onset of instability in 
a weak field is not fullfied. This occurs if the magnetic field is 
very strong or the wavelength of perturbations is small such that 
ujao > to. The latter inquality is equivalent to 



A < A c = 2TTB/VLy 



1.8 x 10' 



cm, (35) 



where A = 27r/fc is the wavelength, B 13 = B/10 13 G, to 3 = 
Sl/10 3 s _1 , pi^ = p/10 14 g/cm 3 . For perturbations, satisfying 
condition (35), the instability can occur if il s < but only in 
the neighbourhood of the cone (or cones depending on the field 
topology) in which the radial component of B is vanishing. The 
opening angle of the unstable region around this cone is of the 
order of SI/ojo (see Eq.(34)). Note that unstable motion in this 
case are perpendicular to the radius as well since gravity sup- 
presses radial motions. The reason why the instability is not en- 
tirely suppressed even by a very strong magnetic field is qual- 
itatively very simple. Perturbations with k = (k r ,0, 0) cannot 
be suppressed by gravity, and such perturbations do not feel the 
stabilizing influence of the magnetic field in the region where 
B r rts as well, because this influence is proportional to the 
Alfven frequency (oc k ■ B). 

It appears that the importance of MRI in core collapse can 
be overappreciated since its growth rate is relatively low and 
reaches the value <~ SI only for perturbations with a wavelength 
<~ A c (see Eq.(35)). Apart from this, gravity strongly suppresses 
development of any perturbations with non-radial wavevectors. 
Only perturbations with k w k r can be unstable but hydrody- 
namic motion for such perturbations has a small radial compo- 
nent and turbulent transport should be inefficient radially. Also, 
it is possible that MRI occurs only in not very extended regions 
of the proto-neutron star that may diminish substantially its ef- 
fect on core collapse. Note that since the wavelength A c is small, 
one needs a very high resolution in numerical simulations to see 
the most rapidly growing modes (Obergaulinger et al. 2009). 
Perturbations with longer A grow substantially slower. 

Our conclusion is in contrast to the results obtained by 
Masada et al. (2007) who considered axisymmetric and nonax- 
isymmetric magnetorotational instability of PNSs taking into ac- 
count dissipative processes. These authors also used a local ap- 
proximation in the stability analysis. In the local analysis, how- 
ever, the shape of perturbations is assumed to be unchanged 
and, therefore, this approach applies only if the growth rate 
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of instability is greater than the rate with which perturbations 
change their shape. In the case of MRI, the growth rate is smaller 
than (or, at maximum, comparable to) s\ VO| whereas perturba- 
tions change their shape because of differential rotation with a 
rate ~ s|Vfi| (see, e.g., Balbus & Hawley 1992 for details). 
Therefore, the results obtained for a nonaxisymmetric MRI in 
the local approximation raise some doubts. As far as axisym- 
metric instability is concerned, Masada et al. (2007) were con- 
fused when identifying different modes in the dispersion equa- 
tion. In the dissipative case, they obtained a dispersion equa- 
tion of the seventh order which describes seven different modes. 
Certainly, different modes can be unstable in different condi- 
tions, but the authors considered only one criterion (Eq.(32) of 
their paper). Unfortunately, this criterion does not correspond to 
MRI but describes a secular instability that is a magnetic analogy 
of the well known Goldreich-Schubert-Fricke (GSF) instability. 
The magnetic analogy of this instability was first considered by 
Urpin (2006) for the case of ordinary stars, and Masada et al. 
(2007) obtained the criterion of the same instability modified 
for the conditions of PNSs. In contrast to MRI, this instability 
is dissipative and disappears if diffusive coefficients go to zero. 
Stratification has a weak impact on this dissipative instability, 
and it can arise in PNSs. Note, however, that diffisive processes 
can reduce the stabilizing influence of stratification on the MRI 
as well. 

As was noted by Masada et al. (2007), diffusion of heat and 
leptons can influence buoyancy if the characteristic diffusion 
timescale is shorter than the buoyant frequency o>o- Since the dif- 
fusion timescale is determined mainly by lepton diffusion which 
is slower than thermal diffusion, this condition is equivalent to 
£fc 2 > ujq where £ is the coefficient of lepton diffusivity. In deep 
layers with a density <~ 10 14 g/cm 3 , the value of £ is ~ 3 x 10 7 
cm 2 s _1 (see Eq.(45) by Masada et al. (2007)). Hence, diffusion 
begins to suppress the stabilizing effect of stratification when the 
wavelength of perturbations is shorter than <~ <~ 10 3 

cm (we suppose uj ~ 10 3 s _1 ). The effects considered in this 
paper are important for longer wavelengths for the standard pul- 
sar magnetic field. 

As was mentioned above, MRI is not the only instability 
caused by differential rotation in proto-neutron stars. The anal- 
ogy of the GSF dissipative instability can occur in both magnetic 
and non-magnetic PNSs. In non-magnetic stars, this instability 
arises if the angular velocity depends on the vertical coordinate 
z (Urpin 2007), and the shellular rotation is a particular case of 
such rotation. Since neutrino transport in PNSs is much more ef- 
ficient than radiative transport in ordinary stars, the dissipative 
instability can be rather fast. If Q z = dVl/dz ~ Q/s, the condi- 
tion of the instability reads 



For perturbations with such wavevector, the growth rate is of the 
order of fi(K T /V)(f2/w ) 2 and is always higher than the growth 
rate of MRI if condition (36) is satisfied. Therefore, dissipative 
instability seems to be more efficient in proto-neutron stars than 
MRI. 
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\sQ z 



fl > 2uj 



(36) 



where v and kt are the coefficients of viscosity and thermal 
diffusivity. The difference to the standard criterion of GSF in- 
stability (<9£1 jdz ^ 0) is caused by the more complex charac- 
ter of neutrino trancport that involves diffusion of both heat and 
the lepton fraction. This instability can arise soon after collapse 
even in dense central regions of the star where convection does 
not occur. It turns out that in all cases, when dissipative insta- 
bility arises, its growth rate is higher than that of MRI. Indeed, 
the maximum growth rate of MRI is <~ SI and can be reached 
for perturbations with wavelength ~ A c . The growth rate of the 
dissipative instability is approximately Kxk 2 ! (0 / 'ojq) 2 . The max- 
imum wavevector for which the instability occurs is ~ yj^ljv. 



